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An efficient computer algorithm is described for the perspective drawing 
of a wide class of surfaces. The class includes surfaces corresponding to 
single-valued, continuous functions which are defined over rectangular 
domains. The algorithm automatically computes and eliminates "hidde?i 
lines." The number of computations in the algorithm grows linearly with 
the number of sample points on the surface to be drawn. An analysis of 
the algorithm is presented, and extensions to certain multi-valued functions 
are indicated. The algorithm is implemented and tested on two different 
computers: a large central computer with hard-copy capability, and a 
small laboratory computer affording interactive use. Running times are 
found to be exceedingly efficient on both machines. Interactive implementa- 
tion of the algorithm, with on-line scope display and view-point control, 
enables effective and rapid examination of a surface from many per- 
spectives. 

I. INTRODUCTION 

The general problem of efficiently and meaningfully displaying 
three-dimensional objects in two dimensions is central to computer 
graphics. In particular, the "hidden line problem" of drawing objects 
in perspective while eliminating line segments not visible from the 
viewing point, is one of long standing. In recent years various algorithms 
dealing with this problem have appeared in literature. 1-7 Some of these 
algorithms entail prohibitively long computation time (relative to 
the number of data points) or have large storage requirements; some 
have both of these undesirable characteristics. 

In this paper a detailed description will be given of a highly efficient 
algorithm for the perspective drawing of an arbitrary surface which 
corresponds to a single-valued continuous function defined over a 
rectangular domain, with elimination of hidden lines. The vantage 
point from which the surface is viewed may be any point not on the 
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surface. Some generalizations of the algorithm to nonrectangular 
domains and to multiple-valued functions will be indicated. 

This algorithm has been implemented on both a large and a small 
computer (Honeywell 6070, DDP 516) with considerably shorter 
computation time than that based upon previous algorithms for drawing 
similar surfaces. A very useful feature of the implementation on the 
DDP 516 (a 16K laboratory machine with 16-bit word length) is that 
the program may be run interactively, allowing the user to select 
varying vantage points and rapidly display different views of a given 
surface. 

II. STATEMENT OF THE PROBLEM 

Given a three-dimensional surface S denned over a rectangular 
domain R, and a vantage point V not on S (Fig. 1), we may assume 
that R is centered about the origin of the x, y-plane and is oriented 
so that its sides are parallel to the z-axis and y-axis, respectively. 
(This may be accomplished without loss of generality by a suitable 
translation and rotation of R, together with V). The plane P, containing 
the perspective image S' of S, is chosen to be perpendicular to the 
line joining V and the origin. A rectangular coordinate system with 
x', y'-axes is chosen on P which preserves the original vertical direction 



V=(v,,v 2 ,v 3 ) 




Fig. 1 — Perspective projection of .surface S on plane P. 
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Fig. 2 — Surface patch of one grid subrectangle. 



of S. That is, the f/'-axis on P must be the projection of the original 
2-axis. This latter choice is necessary, as we shall see later. (Relative 
to the rectangular system on P, the perspective image on P of each 
point of S has a well-defined pair of coordinates. The evaluation of 
these coordinates will be given in the Appendix.) 

The surface S will first be quantized to an M X N rectangular grid 
of points on the domain R. Figure 2 shows the grid lines on R, 
partitioning R into subrectangles. The part of S denned over one 
subrectangle is shown. It will be referred to as a "surface patch" of 
S. After quantization, only the four points of the patch defined at the 
four vertices of the subrectangle are known, and linear interpolation 
of the function will be assumed between adjacent grid points, that is, 
between adjacent points in the same row and adjacent points in the 
same column of the grid. (Thus, S need not be explicitly denned by 
a mathematical function; the data for S may be given by a rectangular 
array of points corresponding to a rectangular grid of points on the 
domain.) The behavior of S in the interior of each subrectangle of the 
grid will be ignored. Thus, each surface patch of S will be represented 
by its four linearly interpolated edges in three-dimensional space. The 
image of these edges on the projection plane is a four-sided polygon 
(Fig. 3). For each surface patch of S, only the visible line segments 
of its edges will be drawn. 

Intuitively, the surface should be thought of as an opaque elastic 
membrane stretched over a rigid frame consisting of all the linearly 
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Fig. 3 — Projection of linearly interpolated surface patch. 

interpolated edges of the surface patches. Then the problem is to give 
a perspective drawing of the frame which eliminates line segments 
of the frame hidden by the opaque membrane from the vantage point. 
(Note that the vertical projection of the frame onto the x, y-plane 
coincides with the grid lines on R.) 



III. BASIC IDEAS OF THE ALGORITHM 

The algorithm for determining the visibility of any given edge is 
based upon two ideas. The first one consists of choosing a particular 
ordering of the surface patches of S so that no part of any 
patch occurring earlier in the ordering is obscured from the vantage 
point by any part of a patch occurring later. The surface is to be drawn 
according to such an ordering, one patch at a time. Then, a line segment 
on any edge of a patch under consideration is hidden from the vantage 
point if and only if its image on the projection plane lies inside a region 
of the plane already covered by the images of earlier patches. 

This particular ordering of the surface patches has the further 
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property that, for a surface S which corresponds to a single-valued 
continuous function over a rectangular domain, the region of the 
projection plane being covered by the emerging image of S will grow 
contiguously as each patch in the ordering is drawn. And because of 
the preservation of the vertical direction of the surface, at each stage 
of the drawing of S, this "hidden" region of the plane can be precisely 
bounded between two piecewise linear continuous functions. The 
determination of the visibility of any edge of a surface patch is thus 
reduced to the problem of deciding which part of its image on the 
projection plane lies between these two functions. This delineation of 
the hidden region of the projection plane, after each patch is drawn, 
by means of two piecewise linear continuous functions, constitutes 
the second idea behind the algorithm. 

Remark: If the y'-axis on the plane of projection is not chosen to be 
the perspective image of the original z-axis, then the image of the surface 
may be "tilted" so that it will be impossible to delineate its boundary 
by using only two functions. 

rV. ORDERING OF THE SURFACE PATCHES 

The fact that it is possible to order the surface patches of S so that 
earlier ones are not obscured by later ones depends on the following 
observation: 

Suppose A and B are any two points of the surface such that A 
obscures B from V. Geometrically, this means that A, B, and V lie 
in a straight line, and the distance from A to V is less than the distance 
from B to V. Then it can be easily shown that the vertical projections 
A , B n , V of the points A, B, V, respectively, onto the x, y-plane 
satisfy this relationship also. That is, they are collinear, and A is 
closer than B to V ■ 

This observation is the key to the ordering of the surface patches. 
First, consider the case in which the vertical projection V = (ui , v 2 , 0) 
of the vantage point onto the x, ?/-plane is southwest of the domain, 
as shown in Fig. 4. [That is, v x ^ x A and v 2 ^ y x , where (x x , y x , 0) 
is the lower left-hand corner of the domain.] An immediate consequence 
of the above observation is that if A = (a x , a 2 , a 3 ) obscures B = 
0>i , b 2 , b 3 ) from V, then a x < 6, and a 2 < b 2 ; i.e., A must be southwest 
of B (given the southwest location of V ). This suggests the following 
rule for determining the relative order of points of £, given the south- 
west location of V : Let A and B be any two points of S. If one of 
them, say A, is southwest of the other (more precisely, if A is south- 
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V=(v l ,v 2 ,v 3 ) 




V =(v,,v 2l o) 



Fig. 4 — Points of the surface with A obscuring B. 



west of B ), then A must precede B. Otherwise, A and B may be ordered 
independently of each other. 

There are various ways to order the surface patches of S so that this 
rule for the relative order of points is satisfied, and hence, so that 
earlier patches are not obscured by later patches. Figure 5 illustrates 
one possible ordering of the subrectangles of R so that the induced 
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Fig. 5 — One possible ordering of subrectangles. 
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ordering of the patches of 5 has this property; i.e., they are ordered 
by rows, from left to right in each row, beginning with the front row 
(relative to the southwest location of V ). It is obvious that, under 
this ordering of the patches, no point of an earlier patch will be obscured 
by a point of a later patch, since no point of a later patch is southwest 
of any point of an earlier patch. 

For other "corner" locations of V (i.e., northwest, northeast, and 
southeast of the domain) analogous orderings of the surface patches 
may be made. Instead of discussing suitable orderings of the surface 
patches for other locations of V , we shall proceed to describe the 
algorithm in detail for the special case in which V a is southwest of the 
domain, and then show how the other cases may be reduced to cases 
involving only "corner" locations of V ■ 

V. THE ALGORITHM IN DETAIL 

Let S be a surface defined over a rectangular domain R of the x, 
?/-plane. Assume S is to be viewed in perspective from a vantage point 
V whose vertical projection V onto the x, y-plane is southwest of R. 
Given an M X N rectangular grid of points on R, for J = 1,2, • • • , M 
and J = 1, 2, • • • , JV, let A(I, J) denote the "sample point" on S 
defined at the grid point on R located at the 7th column and «/th row 
of the grid. For J = 2, • •• , M and J = 2, • • • , JV, let S(I, J) denote 
the surface patch with vertices A (I, J), A (I, J — 1), A (I — 1, «/), 
A(I — 1, J — 1). (See Fig. 2.) By above considerations, the patches 
may be ordered as follows, and drawn according to this ordering, 
one at a time: 

S(2, 2), 5(3, 2), • • • , S(M, 2), 

S(2, 3), 5(3, 3), • • ■ , S(M, 3), 



S(2,N),S(3,N), ■■■ ,S(M,N). 

Now, consider the sequence of points A(l, J), J = 1, ■ ■ • , N, defined 
along the first column of grid points. The line segments joining these 
points in consecutive order will be called the "leading left edge" of S. 
Similarly, the line segments joining the sequence of points A (I, 1), 
I = 1, • • • , M in consecutive order will be called the "leading front 
edge" of S. As a consequence of the observation made at the beginning 
of Section IV, both leading edges are entirely visible from V. Thus, 
we may begin by drawing all the line segments of the leading edges. 
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Then the drawing of each successive patch S(I, J) in the ordering may 
be completed by adding the visible segments of just two of its four 
edges, namely its "back edge," joining A(I, J) and A(I - 1, J), and 
its "right edge," joining A(I, J) and A(I, J - 1), since the other two 
edges will have already been considered in connection with earlier 
patches or with one of the leading edges. 

Figure 6a gives an example of a surface denned over a 5 X 3 grid 
of points. The 5X3 array of sample points on the surface gives rise 
to two rows of patches, four in the front row and four in the back row. 
Figure 6b shows the conjunction of the leading left edge and the leading 
front edge of this same surface. 

Figures 7a through 7h show successive patches being added to the 
drawing of the surface. With each patch that is being added, L, denotes 
its "back" edge, and L 2 denotes its "right" edge. Dotted lines indicate 
hidden segments which are not drawn. 

As indicated earlier, the determination of the visibility of each new 
edge under consideration will involve two piecewise linear continuous 
functions. Intuitively, one of these functions, Max, will be used to 
define the path of the upper perimeter of the region of the projection 
plane thus far covered, and the other function, Min, will be used to 
define the path of the lower perimeter. At this point, it would be useful 
to identify the grid of the plotting device with a rectangle in the pro- 
jection plane which circumscribes the perspective image of the surface. 
(This corresponds to a scaling and translation of the coordinates of 
the image points so that they fall within the dimensions of the plotting 
grid.) Then, in order to optimize the precision in delineating the hidden 
region of the projection plane, or equivalently, the hidden region of 
the plotting grid, the number of breakpoints of each of the functions 
Max and Min is chosen to coincide with the number of horizontal 
units of the plotting grid. 

Before plotting begins, we artifically set Max (fc) = Min (k) = — 1, 

ACl.3) 

A(3,1) 





A(5,1) 
A(4,1) 



(a) (b) 

Fig. 6 — Surface and its leading left and front edges. 
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A(2,2) -L 




A(3,2) ^|_ 2 
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(d) ^ (h) 

Fig. 7 — Successive patches of surface being drawn. 



for k = 1, 2, • • • , K x , where K z is the horizontal dimension of the 
plotting grid. Then, after each line segment is drawn, Max or Min 
(possibly both) will be modified to reflect the change in the shape of 
the emerging image of S. We begin with the line segments on the leading 
edges of *S. Let A and B be any two consecutive "sample points" of 
the leading left edge or of the leading front edge, that is, any two 
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consecutive points in the following sequence (see Fig. 2): 

A{\, N), A{\, N -1), ■■■ , A(l, 1), A(2, 1), • • • , A(M, 1). 

Let A x , A„ be the (scaled and translated) plotting coordinates of A, 
and B z , B u the plotting coordinates of B. First, the line segment joining 
the points (A x , A u ) and (B z , B v ) is drawn. Then both Max and Min 
are redefined between A z and B z by setting 

Max (A z ) = Min (A T ) = A u 

Max (B x ) = Min (B x ) = B u 

and linearly interpolating both Max and Min between A x and B x . 
(In Fig. 8, K x and K u refer to the dimensions of the plotting grid, in 
the horizontal and vertical directions, respectively.) 

After the above procedure is carried out for every pair of consecutive 
points in the sequence, the leading edges will be completely drawn, 
and their path will be defined by Max and Min between the two end- 
points. Note that the region bounded between Max and Min has area 
zero, reflecting the fact that no patch has yet been drawn. 

We now proceed to draw the surface patches. Each S(I, J) in the 
ordering, beginning with S(2, 2) will be uniformly processed as follows: 
Let L y denote the edge between A (I, J) and A (I - 1, J), L 2 denote 
the edge between A(I, J) and A(I, J - 1). For each of these edges, 
we must first determine its visible segments, then draw only those 
segments, and modify Max or Min to reflect the change in the shape 
of the hidden region (the entire edge may be hidden, of course). 

First, let L = L, and A = A(I, J), B = A(I - I, J), the endpoints 
of Li . The following criterion is used to determine the visibility of 
any point C on L: If C z , C„ are the plotting coordinates of C, then 
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Fig. 8 — Linear interpolation of Max and Min. 
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C is visible if and only if C„ ^ Max (C x ) or C u ^ Min (C z ). This visibility 
criterion is applied to each of the endpoints A and B. (Let A z , A v 
be the plotting coordinates of A; B x , B„ those of B.) There are four 
possibilities: 

(i) A, B are both visible. Then all of L is assumed to be visible. 
(This assumption is made to reduce computation time. It may happen 
that some segments of hidden lines will be incorrectly drawn if the 
surface has sharp "spikes" in the foreground. This problem can be 
circumvented by starting with a finer grid on R.) The line segment 
joining (A x , A u ) and (B z , B„) is drawn. If A„ ^ Max (A z ), Max is 
linearly interpolated between A z and B z by setting Max (A z ) = A u , 
Max (B z ) = B u . Otherwise, we must have A u ^ Min (A z ), and Min is 
similarly modified between A z and B z , with Min (A z ) = A„ , Min (B z ) = 
B„. 

(u) A and B are both hidden. Then for the same reasons as in case 
(i), we assume all of L is hidden and no line segment is drawn. Max 
and Min are left unchanged, and we proceed to the next edge. 

(iii) A is hidden and B is visible. A search is made along the line 
joining A and B, starting at A, for the first visible point C (discrete 
steps are taken in the horizontal direction of the plotting grid). The 
segment joining C to B is assumed to be entirely visible and is drawn. 
If C u ^ Max (C z ), Max is linearly modified between C z and B z with 
Max (C z ) = C u , Max (B z ) = B u . Otherwise, Min is modified between 
C x and B z . 

(iv) A is visible and B is hidden. A search is made along the line 
joining A and B for the first visible point C, starting from the hidden 
point B. The rest is analogous to case {Hi). 

For each surface patch S(I, J), this procedure is carried out with 
L = Li and B = A(I — 1, J); then the procedure is repeated with 
L = L 2 and B = A (I, J — 1). This completes the description of the 
algorithm for the special case in which V , the vertical projection of 
the vantage point onto the x, i/-plane, is southwest of the domain R. 

If V is in any of the other "corner" regions of the x, y-plane (indicated 
by NW, NE, and SE in Fig. 9a), we may rotate the data defining the 
surface (either the mathematical function or the rectangular array 
of sample points) by 90, 180, or 270 degrees, respectively, and apply 
the algorithm as described for V in the SW region of the x, ?/-plane. 
Another approach is to change the ordering of the surface patches and 
to define new leading edges to suit the other locations of V . 

If V Q is in any of the regions directly west, north, east, or south of R 
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(a) (b) 

Fig. 9 — Locations of vertical projection of vantage point. 

(Fig. 9b), we may partition R into two subrectangles Ri and R 2 
(Fig. 10a), and apply the algorithm to each of B\ , R% , since V is in 
a "corner" location relative to each subrectangle. The two parts of 
the surface, corresponding to Ri and R 2 , will have to be "pieced to- 
gether" along the line of partition, and this may be done without 
"flaws" by augmenting the sample points of the surface with the set 
of points denned at the intersection of the partition line with the grid 
lines on R. 

Finally, if V lies inside the domain R, we may partition R into four 
subrectangles (Fig. 10b) and apply the algorithm to each of R x , R» , 
R 3 , R* , with V in a "corner" location. 

Remark: Note that under the assumption that the plane of projection 
P be perpendicular to the line joining V and the origin, if V coincides 
with the origin, then the projection of the z-axis from V onto P is a 
point. In this case, the y'-axis on P must be chosen to be the image 
of some other line. See the Appendix. 

Figures 11, 12, and 13 are sample drawings generated on the Honey- 
well 6070 computer with the Stromberg-Carlson 4060 microfilm plotter. 
Each surface was evaluated at a 64 X 64 grid of points, for a total of 
4096 data points. The average run time was 9 seconds for each surface. 
(Part of the 9 seconds was used in computing the function values and 
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Fig. 10 — Partitioning of domain for non-corner locations of V . 
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plotting coordinates at each point of the 64 X 64 grid. In fact, these 
computations were performed twice, once for determining the scale 
factors and again for the actual plotting of the surface.) Total storage 
for generating microfilm output of each surface was 16K. 

It is clear that using this algorithm, computation time should vary 
linearly with the number of data points, since each data point is pro- 
cessed by evaluating its plotting coordinates IX, IY, and then com- 
paring IY with just two values, Max (IX) and Min (IX). 

Economically generated movies have also been produced on both 
the Honeywell 6070 and the DDP 516 computers. These movies show 
surfaces in rotation and surfaces undergoing continuous change in 
shape. The average run time on the Honeywell 6070 for making such a 
movie with a 32 X 32 grid and viewed from 240 distinct vantage points 
was 7.4 minutes, or about 1.85 seconds per frame. The corresponding 
run time on the DDP 516 was 7.5 seconds per frame. (For these movies, 
the scale factors were precomputed so that function values and plotting 
coordinates were evaluated only once for each frame of the movie.) 

For a rough comparison of computation times for hidden-line elimina- 
tion drawings based upon other algorithms (and implemented on other 




Fig. ll—/(i, y) = (sin (x) - 1) (sin (y) - 1), -» g x,y g w. 
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Fig. 12— f(x, y) = sin (x + y)/(\ + (x - y) 2 ), -2ir £x,y ^ 2tt. 



computers) see Refs. 1 through 7. In particular, see Refs. 2 and 3 for 
a comparison with two other algorithms which deal with single-valued 
functions denned over a rectangular domain. Using the algorithm 
described in Ref. 2, computation time increases quadratically with 
the number of data points, since the visibility of each data point is 
determined by an exhaustive comparison with every other data point. 
Although computation time increases linearly using the algorithm 
of Ref. 3, the storage required is large in order to obtain high-quality 
drawings. This is because the "plotting page" is subdivided into rec- 
tangles, and, of all possible line segments of the surface whose images 
lie in any given rectangle, only the one closest to the vantage point is 
drawn in that rectangle. Thus, in order to avoid a "sketchy" drawing, 
a fine subdivision of the "plotting page" is necessary and a large array 
is required to store the information in each rectangle of the subdivision. 
(The algorithm of Ref. 3, however, is applicable to a larger class of 
surfaces than the algorithm described in this paper. For example, it 
can be used to draw surfaces corresponding to functions which are 
not single-valued, e.g., intersecting cylinders.) 
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VI. GENERALIZATIONS 



6.1 Convex Domains 



The algorithm may be generalized to an arbitrary surface S defined 
by a single- valued continuous function over a convex domain R, by 
artificially extending the denning function / to a rectangular domain 
R' containing R. The ideas of the algorithm are then applied to the 
surface S' defined by the extended function /' over R', with the following 
modifications: First, the sample points along the leading edges of S' 
must be replaced by a subset of points defined at the intersection of 
the boundary of R with the grid lines on R', and special consideration 
must be given to drawing the lines connecting these points. Second, 
a test must be made of sample points of S' to determine if they belong 
to S, so that only line segments that are part of S are drawn. One way 
to accomplish this is by setting f(x, y) = z for all (x, y) in R' — R, 
where z Q is a number outside the range of values of / on R. (Convexity 
of the domain is necessary in order to bound the emerging images at 
each stage of the drawing between two functions.) 
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Fig. 13— f(x, y) = cos (i 2 + 7/i)/e - 2 <*'-**'> -2.5 g x, y g 2.5. 
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6.2 Multi-Valued Functions 

The algorithm may also be extended to certain multi-valued functions. 
For example, if the surface is a sphere, it is first partitioned by a hori- 
zontal plane into two single-valued pieces, and if the vantage point 
is above the horizontal plane, the algorithm is first applied to the upper 
hemisphere, and then, without reinitializing Max or Min, the algorithm 
is applied to the lower hemisphere. 

6.3 Bivariate Histograms 

With minor modifications the algorithm may be adapted to the 
perspective drawing of bivariate histograms, with elimination of hidden 
lines. Figures 14a through 14d show four different perspectives of a 
bivariate histogram generated by the same data. The computation 
time on the Honeywell 6070 was 1.6 seconds for each of the four drawings. 




Fig. 14 — Bivariate histograms: (a) V in SW corner, (b) V in SE corner (c) V . 
in NE corner, (d) V in NW corner. 
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APPENDIX 



Evaluation of Rectangular Coordinates of Image Points 

Let S be a surface defined over a domain centered at the origin of 
the /, //-plane, and let V = (u, , v 2 , v a ) be a point not on S. In order 
to simplify the calculation of the coordinates of image points and hence 
reduce computation time, a special plane P*, normal to Y and passing 
through the origin, will be used as the plane of projection. (In Fig. 1, 
a plane P parallel to P* is shown. Note that the image of S on P differs 
from the image of S on P* only by a scale factor.) 

Let A = (ai , a 2 , a 3 ) be any point on S, and let B = (by , b 2 , b 3 ) 
denote its image point on P*. The plotting coordinates of A depend on 
the two-dimensional c oordinates (u , v) of B with respect t o an orthogonal 
system on P*. If X = (.t, , x 2 , x 3 ) and Y = (y x , y 2 , ?/ 3 ) are the positive 
unit vectors of the coordinate system on P*, then the scalars u and v 
must satisfy the following relationship: 

uX + vY = B. (1) 

Since B is collinear with the points A and V, then 

B = V + r(A — Y), for some scalar r. (2) 

Thus, eq. (1) may be replaced by 

uX + vY = V + rA - rV. (3) 

Taking the dot product of both sides of eq. (3) with X, we have 

uX-X + v?-X = YX + rAX - rYX. 

Since XX = 1 and YX = V X = (V-X = because X on P* 
and P* normal to V), we have 

u = rA-X. (4) 

Similarly, taking the dot product of both sides of eq. (3) with Y, we 
have 

v = rAY. (->) 

The scalar r is determined by taking the dot product of both sides of 
eq. (2) with V. Since B lies on P*, then BY = 0, and we have 

r = V-V/V-(V - A). (6) 

It remains to determine the coordinates of the unit vectors X, Y. 
As mentioned in Section II, the ?/-axis on P* must be the perspective 
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image of the original z-axis. Note, however, that if V is on the z-axis, 
then the perspective image from V of the z-axis is a single point, which 
of course cannot be used as an axis line. In this case, we must com- 
promise by either moving the vantage point slightly off the z-axis 
or choosing the y'-axis on P* to be the image of a slightly "askew" 
z-axis. For the rest of this discussion we assume that V is not on the 
z-axis. 

Let Z = (0, 0, 1). Then V, Y, Z lie on the same plane, because Y 
is the image from V of some scalar multiple of Z. Thus Z is a linear 
combination of V and Y. Since XV = and X- Y = 0, then XZ = 0. 
The normalization of any vector X' satisfying the equations X'V = 
and X'Z = will be a unit vector on the s'-axis of P*. X' = 
( — v 2 , Vi, 0) satisfies these two equations. Let d, = -\/v\ + v\ . Then 

X = (-!>,/<*! , v 1 /d 1 , 0) 

is the positive unit vector on the rc'-axis of P*. The coordinates of Y 
are similarl y determined from the equations Y ■ X = and Y-V = 0. 

Let d = Vvl + v\ + vl . Then 



Y = (-V&3 , —v 2 v 3 , di)/dd t 



is the positive unit vector on the y'-axis of P*. Substituting these values 
of X and Y back into eqs. (4) and (5), we have 

u = r(a 2 Vi — a 1 v 2 )/d l (7) 

v = r{—a x v x v 3 — a 2 v 2 v 3 + o 8 d?)/(ddi), (8) 

where r = d 2 /(d 2 — (o^ + a 2 v 2 + a 3 v 3 )), as denned in eq. (6). 

The above expression for v is algebraically equivalent to the following 
one, which involves half as many arithmetic operations: 

v = (v 3 + r(a 3 — v 3 ))-d/di. (9) 

Thus, given a point A on S, the rectangular coordinates (w, v) of 
its image point are evaluated according to eqs. (7) and (9). 
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